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Angular multiscale statistics of Lagrangian trajectories in turbulence 
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The angle between subsequent particle displacement increments is evaluated as a function of the 
timelag in isotropic turbulence. It is shown that the evolution of this angle contains two well-defined 
power-laws, reflecting the multi-scale dynamics of high-Reynolds number turbulence. The proba¬ 
bility density function of the directional change is shown to be self-similar and well approximated 
by an analytically derived model assuming Gaussianity and independence of the velocity and the 
Lagrangian acceleration. 
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Advances in experimental devices and numerical sim¬ 
ulations over the last two decades have opened the way 
to a Lagrangian characterization of turbulent flows [ij 
0- The structural description of the statistical dynam¬ 
ics of turbulence has thereby shifted from the investi¬ 
gation of spatial correlations of instantaneous velocity 
fields to the study of temporal correlations along fluid 
particle trajectories. In the Lagrangian reference frame, 
the spatio-temporal complexity of turbulence manifests 
itself through the spiraling chaotic motion of fluid par¬ 
ticles, changing direction at every timescale. This di¬ 
rectional change of Lagrangian tracers, as a function of 
the timelag between two observations, is the subject of 
the present work. Instantaneous measures of the curva¬ 
ture in turbulence have been investigated in the past five 
years for academic turbulent flows, both in three m 
and in two space dimensions 0,0. Curvature is dom¬ 
inated by the small-scale structures and contains only 
little information on the multiscale dynamics of turbu¬ 
lent flows. Multi-scale dynamics can be measured by 
Lagrangian structure functions [ 3 , 0 , but those do not 
contain any direct information on the curvature of the 
trajectories. 

The related measure which represents the coarse 
grained curvature over a time interval was only recently 
introduced by Burov et al. 0. More precisely, in this last 
work the directional change of a particle was introduced, 
and the characteristics of this new measure were investi¬ 
gated in various types of random walks. In the present 
work, we will show how this measure can characterize 
the time-correlation of the direction of a fluid particle 
in a turbulent flow. In particular will we show how the 
multi-scale character of a turbulent flow can be revealed 
by considering the timelag dependence of the directional 
change. 

We define the Lagrangian spatial increment as 

SX(x 0 ,t, t ) = X(x 0 ,t) - X(x 0 ,t - t ) (1) 
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FIG. 1: Top: definition of the angle between subsequent La¬ 
grangian particle increments. Bottom: short time evolution 
and definition of the lengthscales l± and l \\. 


increments, introduced in js], is 


COsiOlf Til = SX ( X 0^M)-5X{x 0 ,t + T,T) 

\SX(x 0 ,t,T)\ \SX(x 0 ,t + t , t )\ 


( 2 ) 


The angle is illustrated in Figure [T] (top). Rather than 
considering its instantaneous evolution, its averaged ab¬ 
solute value is of particular interest in an isotropic ran¬ 
dom velocity field. The ensemble average will be denoted 
in the following by 


where X(x$ ,t) is the position of a fluid particle at time 
t , passing through point Xq at the reference time t = to 
and advected by a velocity field u, i.e. dX jdt = u . The 
cosine of the angle 0(<, r) between subsequent particle 


0(r) = <|e(i,r)|>. (3) 

We omitted the time-dependence since we will consider 
statistically stationary flow in the following. For short 
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FIG. 2: The average angle 9 as a function of the timelag r 
normalized by the Kolmogorov time-scale tk■ In the inset the 
compensated angle, 6(t) = 6(r)a u / (er ) 1 ' 2 is plotted. 


time lags, 9{t) should be close to zero, whereas for times 
long compared to the correlation time associated with the 
spiraling motion 9(r) should tend to 7 t /2 by symmetry. 

For short times the instantaneous angle 0(r, t) is re¬ 
lated to the curvature k (see Figure [T] (top)) by the rela¬ 
tion 


n(t) = lim 


|0(f,r)| 


r-y 0 2r||M(t)|| ' 


( 4 ) 


with u being the velocity. How the angle varies in be¬ 
tween the short and long-time limits is the main subject 
of the present work and we will show that the depen¬ 
dence of 6{t) on the timelag contains the signature of 
the multi-scale dynamics of a turbulent flow. 

The database used to investigate the behaviour of 9{t) 
is described in l,E3. The simulation was carried out us¬ 
ing standard pseudo-spectral techniques, following 8 . 10 6 
fluid particles in a statistically stationary isotropic turbu¬ 
lent flow during 5.8 integral timescales in a periodic cube 
of dimension 27r. The resolution is 1024 3 gridpoints. The 
integral timescale is 2.1 and the Kolmogorov timescale 
tk = iy/e) 1 ! 2 = 0.036, where e = 0.31 is the mean 
dissipation rate and v = 4.10 ~ 4 the kinematic viscosity. 
The Lagrangian integral timescale is of the order of the 
Eulerian integral timescale. The Taylor-scale Reynolds 
number is R\ = 225. 

Figure [5] shows 9(r) in double-logarithmic representa¬ 
tion. The angle increases monotonuously from zero to 
7t/2, and this latter value is approached for values of r 
of the order of the Lagrangian integral timescale. Two 
power-laws can be identified in this graph, with a cross¬ 
over around twice the Kolmogorov timescale. The origin 
of these power-laws will now be elucidated. 

For our phenomenological explanation, we consider 
high-Reynolds-number isotropic turbulence, containing 
flow structures on a wide range of different scales. We 
consider short timelags r <C T, where T is the Lagrangian 


integral timescale of the flow. In this limit, the angle 
0(t,r) can be approximated using a Taylor-expansion 
by, 

^«|tan(0/2)|«|0/2|. (5) 

Ml 


where l± and 1 1 | are shown in Fig. [l] (bottom), and corre¬ 
spond to the absolute values of the distance travelled par¬ 
allel with, and perpendicular to, the initial displacement 
increment, respectively, over a time interval 2r. The val¬ 
ues of 2 1 | and l± can be estimated, again using a Taylor 
expansion, to be 

Z|| « 2[/(f, t)t l± ~ 2r 2 aj_(f, t), (6) 

with U(t,r) and a±(t,r) the absolute values of the ve¬ 
locity and the acceleration perpendicular to the velocity, 
respectively, coarse-grained over a time r along the fluid 
particle trajectory. Without loss of generality, we will 
write t/(f,r) and a±(t, r) as 


U(t,r) = a u (T)^ u (t,T) a_i_(t,r) = cr a (r)i; a (t, r), (7) 


where ct 2 (t) and ct 2 (t) are the variance of the coarse¬ 
grained velocity and perpendicular acceleration, respec¬ 
tively. The quantities £ u (t, r) and ^ a (t,r) are positive 
random variables with unit mean value and unit mean 
variance. We thereby obtain, 


|0(t,r)| 


0 _ °a (r)g a (t,r) 
<T u (r)^ u {t,T)' 


( 8 ) 


We assume the velocity and the acceleration independent, 
a reasonable assumption at very high Reynolds numbers, 
as long asrCT. Without coarse-graining this assump¬ 
tion was also used in reference @ to model the curvature 
in isotropic turbulence. Using this assumption we find 


9{t) 


2 r 


<Tg( T ) 

<Tu(t) ' 


( 9 ) 


Since for t T the velocity is roughly constant over 
the time-interval, <j u (t) w f7 rms . However, aj_ is dom¬ 
inantly determined by the small scales and fluctuates 
rapidly. Only for r small with respect to the smallest La¬ 
grangian time-scale, the Kolmogorov scale, can we con¬ 
sider tr a (r) ~ (fli)rms 5 be., independent of r. For these 
very short time-scales, we have thus 

9(t) k, 2 t— for t <C tk- ( 10 ) 

where a a and a u are the total rms perpendicular accel¬ 
eration and velocity, respectively. The linear relation be¬ 
tween 9{t) and r is well observed in Figure [2] We can 
further express this in terms of quantities which are easy 
to determine experimentally. Assuming classical scaling 
a the acceleration variance is given by the relation 


a 


2 

a 


e 3/2 

^72’ 


( 11 ) 
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a relation which can be refined to take into account inter- 
mittency corrections [l2l - fl4j . Omitting these corrections, 
expressions m and m yield, 

0( T ) ~ 7^ r1 x 2 for T < t k ■ ( 12 ) 

At timescales larger than tk, but smaller than T, i.e., 
in the inertial interval, the above approximations to ob¬ 
tain © are still valid. However, the subsequent ap¬ 
proximation, that <j 0 (t) is independent of r is not valid 
anymore. Indeed, the perpendicular acceleration fluctu¬ 
ates rapidly in time, on a time-scale of the order of tk - 
Coarse-graining the acceleration over an interval r > tk, 
the influence of the more rapidly fluctuating scales is fil¬ 
tered out. Indeed, even if their contribution to the rms 
acceleration is dominant, if the coarse-graining is per¬ 
formed before considering the norm, positive and nega¬ 
tive contributions will cancel each other. The remaining 
variance will be predominantly caused by scales with a 
time-scale larger than, or comparable to r. Following 
classical Kolmogorov phenomenology mill , the accel¬ 
eration induced by inertial range structures with typical 
timescale r will be of the order 

<7a(r) ~ (e/r) 1/2 . (13) 

This estimate is obtained by neglecting the viscous con¬ 
tribution to the acceleration, a reasonable assumption 
even near the dissipation range scales [T 3 J, and realiz¬ 
ing that the acceleration is due to pressure forces, which 
satisfy, at inertial range scales to a good approximation 
Kolmogorov-scaling (Tt], |l8[ . The scale of such eddies is 
proportional to 

l(r) ~ t 3/ V /2 . (14) 

The reciprocal dependence of the acceleration variance 
on r in expression m illustrates that the smallest scales 
are most efficient in accelerating the fluid particles. Af¬ 
ter the influence of the scales smaller than r is removed 
by the coarse graining, the remaining dominant contribu¬ 
tion is caused by the smallest scales still present, i.e., with 
timescale r. It is therefore those scales, with correlation¬ 
time r which will deviate particles from their trajectory 
over a lengthscale of the order of the correlation-length 
of the structures. This phenomenological picture is il¬ 
lustrated in Fig. [I] bottom, where it can be understood 
intuitively that scales of the size I < I(r) are too small 
to efficiently contribute to a perpendicular displacement 
averaged over a time-interval r. 

Combining (fTSl) and © we obtain in the inertial range 

e 1 / 2 / T N 1/2 

9{t) ~ t 1 / -~ ( — ) for t k <C t -C T. (15) 

Again, this scaling is observable in Figurc[2j even though 
the power-law is less well present than in the dissipa¬ 
tion range. This is better appreciated by considering the 
compensated angle, 9(t) = 9(r)a u / (et) 1 / 2 , plotted in 



(a) 



cos(©) 

(b) 


FIG. 3: PDFs of (a) 0 and (b) cos(0) for different timelags. 
Solid black lines indicate the long time asymptotic form of 
the PDFs. 


the inset of the figure. The slow emergence of inertial 
ranges with the Reynolds number in Lagrangian statis¬ 
tics is fairly common El and is was recently even argued 
that they might be non-existent [20]. In the present case, 
the emergence of a plateau is undeniable. This might 
be because the inertial range scaling of the mean-angle 
9{t) is not directly related to the Lagrangian structure 
functions. Indeed, the scaling is induced by considering 
the coarse-grained Lagrangian acceleration, a quantity of 
which the scaling is related to that of the Eulerian pres¬ 
sure gradient. 

The above arguments and results considered the av¬ 
erage value 9{t) only. Further information, in particu¬ 
lar on higher order moments, is contained in the prob¬ 
ability density function (PDF) of the instantaneous an¬ 
gle and its evolution with r. Those PDFs of the angle 
0(f,r) and its cosine are shown in Figure [3] It is ob¬ 
served that the P T (0) for small r consists of a peak near 
zero, whereas for long-times a symmetric distribution be¬ 
tween 0 and tt is obtained. This latter distribution corre¬ 
sponds to the distribution between two randomly chosen 
vectors in three-dimensions. Its distribution is given by 
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where we used expression J5]). If we assume a_L to sat¬ 
isfy a Gaussian distribution, which is only a good ap¬ 
proximation for the core of the PDF, and if we further 
assume u to be uncorrelated with aj_, and also multivari¬ 
ate Gaussian [5j , then both £ 2 and £ 2 satisfy a y-squared 
distribution. For a given velocity vector in 3D, having 3 
degrees of freedom, the perpendicular acceleration is con¬ 
fined to the plane perpendicular to the velocity and is a 2- 
component vector. The ratio of two properly normalized 
independent x 2 -distributed quantities with n, m degrees 
of freedom, respectively, is given by an F n , m Fischer- 
distribution. We expect 1 — cos(0) therefore to be given 
by an F 2 ,3 distribution. More precisely, 

7(' r )-Pi-co 8 (0(t,r))(a;/7(' r )) =F 2 , 3 (x), (17) 



(b) 

FIG. 4: (a) PDFs of 1 — cos(0). (b) The same PDFs normal¬ 
ized and compared to the analytical prediction. 


Poo(0) = sin(0)/2. The distribution of the cosine of the 
angle between two random vectors is thus equidistributed 
so that Poo(cos(0)) = 1/2. It is observed that these two 
long-time distributions are approached for long times in 
Figure [3] We will now show how we can predict the 
short-time, small 0(r, t) behaviour of the PDFs. In par¬ 
ticular will we consider the PDF of 1 — cos(0). Since the 
cosine of small deviations in 0 gives values near unity, 
1 — cos(0) directly measures the magnitude of the di¬ 
rectional change. In addition, it is easy to compare the 
PDF of this quantity to the long-time limit consisting 
of a straight line. The PDFs of 1 — cos(0) are shown 
in Figure |4|a) in double logarithmic representation. To 
explain their shape and their evolution with r, we use a 
Taylor expansion for small 0, 


1 — cos(0(f, t)) 



T ) 2 


2t 2 




(16) 


where y r = 0(r) 2 /3 and 0(r) is shown in Figure [2] It is 
observed in Figurc[4jb) that the agreement with the pre¬ 
diction is fairly satisfactory considering the assumptions 
we made in the derivation of the shape of the PDF. Note 
that no adjustable parameters were used to fit the PDF 
to the P-distribution. 

The results obtained in the present investigation show 
that the time-series of the Lagrangian position can re¬ 
veal the inertial range structure of turbulence through 
the timelag-dependence of the quantity 0{t). In partic¬ 
ular do we show how Kolmogorov’s inertial range theory 
is linked to the angular statistics of Lagrangian fluid par¬ 
ticle trajectories. 

The present framework will allow experimentalists to 
verify the scaling of Lagrangian statistics in very-high- 
Reynolds numbers flows, even if the measurement tech¬ 
niques are not sufficiently rapid to resolve down to the 
Kolmogorov scale. Indeed, no measurements of the in¬ 
stantaneous velocity or acceleration are needed, only La¬ 
grangian position measurements sufficiently sampled to 
resolve the intertial range timescales. A further come- 
out of this investigation are the scale-dependent measures 
for the mean-angle and the probability density functions, 
which will allow to more accurately model the topology 
of Lagrangian trajectories in dispersion models. 

The measure we investigated in the foregoing allows 
a different angle of attack on the simultaneous charac¬ 
terisation of the multi-scale character of turbulence and 
the scale dependent curvature of Lagrangian fluid parti¬ 
cle trajectories. In this light an interesting perspective is 
to clarify the link between the current work and the re¬ 
sults obtained using the recently introduced lon gitu dinal 
and transversal Lagrangian structure functions [21(. 
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